Damped bead on a rotating circular hoop - a bifurcation zoo 
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The evergreen problem of a bead on a rotating hoop shows a multitude of bifurca- 
tions when the bead moves with friction. This motion is studied for different values 
of the damping coefficient and rotational speeds of the hoop. Phase portraits and 
trajectories corresponding to all different modes of motion of the bead are presented. 
They illustrate the rich dynamics associated with this simple system. For some range 
of values of the damping coefficient and rotational speeds of the hoop, linear stability 
analysis of the equilibrium points is inadequate to classify their nature. A technique 
involving transformation of coordinates and order of magnitude arguments is pre- 
sented to examine such cases. This may provide a general framework to investigate 
other complex systems. 
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I. INTRODUCTION 

The motion of a bead on a rotating circular hoopP shows several classes of fixed points 
and bifurcationpJSl. It also exhibits reversibility, symmetry breaking, critical slowing down, 
homoclinic and heteroclinic orbits and trapping regions. It has been shown to provide a 
mechanical analogue of phase transition^. It can also operate as a one-dimensional pon- 
deromotive particle trapP. The rigid pendulum, with many applications, can be considered 
a special case of this systemP^. 

In this article we examine the motion of a damped bead on a rotating circular hoop. 
Damping alters the nature of the fixed points of the system, showing rich nonlinear fea- 
tures. The overdamped case of this modeP^ and a variant involving dry frictiorP^ has been 
previously studied. 

For certain values of the damping coefficient and the rotational speed of the hoop, linear 
stability analysis predicts a line of fixed points and some of the fixed points appear as 
degenerate nodes. However, such fixed points are borderline cases, sensitive to nonlinear 
terms. By transforming to polar coordinates and employing order of magnitude arguments 
we analyze these borderline cases to determine the exact nature of these fixed points. To 
our knowledge, such analytical treatment does not appear in literature. The basic equations 
obtained for this system are quite generic and arise in other systems (e.g. electrical systems) 
as well. Hence, our technique may serve as a framework for investigating other more complex 
nonlinear systems. 



II. THE PHYSICAL SYSTEM 



A bead of mass m, moves on a circular hoop of radius a. The hoop rotates about its 
vertical diameter with a constant angular velocity u. The position of the bead on the hoop 
is given by angle 9, measured from the vertically downward direction {—z axis), and is the 
angular displacement of the hoop from its initial position on the x-axis (Figure [T]). 

The Lagrangian of the system with no damping is, 
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L(9, 6) = ^ (6*^ + sin^ 9) + mga cos 9 , where a; = is a constant. (1) 
Using the Euler-Lagrange equation, the equation of motion is obtained as, 

9 = sin 9{uj'^cos9 - g/a). (2) 



FIG. 1. Schematic figure of a bead sliding on a rotating hoop showing the angles 9 and (p. 



To include friction, a term —b9 is introduced in ([2j) as, 

e = sme{uj\ose - g/a) -W , (3) 

where b is the damping coefficient. We identify u"^ = g/a as the critical speed of rotation of 
the hoop, and write k = w^/w^, fi = b/uc- Defining r = cjc ^, ([3]) may be made dimensionless 
by changing from t to r, 

6" = sm6(kcos6 - 1) - fiO' , where — ■ = 6" . (4) 



For phase plane analysis, we define a new variable 9i = 9 , and write (|4j) as, 

9' = 9, (5) 
^i' = -sin0(l-A;cos^)-/i^i . (6) 

The parameter k can take only positive values whereas /i may be either positive or negative. 

Due to the symmetry of the hoop about its vertical axis, (§ and g remain invariant 
under the transformations 9 — )■ —9 , 9i — t- —9i. This implies that alternate quadrants of the 
9 — 9i plane have similar trajectories. Similarly, it is easily verified that if {9{t),9i{t)) is a 
solution for positive damping (yU > 0), then for negative damping (yU < 0), {9{—t), —9i{—t)) 
and {—9{—t),9i{—t)) are two solutions. The phase portrait of the system for negative 
damping will just be the reflection of the positive damping phase portrait with the arrows 
reversed. Hence we confine our attention to G [0,7r] and /i > 0. 
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When there is no damping^^, the fixed points are at (0,0) and (vr, 0) for < /i: < 1, 
whereas for A; > 1, an additional fixed point appears at {Qi = cos~^(l//c), 0). Damping 
changes the nature of fixed points and not their number or location. 

III. NATURE OF THE FIXED POINT (0, 0) 

The Jacobian matrix at (0, 0) is obtained by Taylor expanding ([s]) and ^ about (0, 0) 
and retaining the linear terms. 



J(0,0)= . (7) 




Let r and A denote the trace and determinant of the above matrix. 

1. When k > 1, both F and A are negative. The fixed point is a saddle with eigenvalues 
and eigenvectors given by, 

/ 



M,2 — 7. ) Vi 2 



1 

\^(-/i±ei)/2 



where = a//!^ + 4(A; — 1). Saddles are robust and do not get perturbed by non- 
linearities. Thus, (0, 0) will remain a saddle even if nonlinear terms are taken into 



account (see Figures M M and 16 ) 



For /i = 0, both Ai and A2 equal \/k — 1. As A; — )■ 1+, Ai — i- 0+ and A2 — > — /i^, which 
means that the saddle will start looking like a line of fixed points along the direction 
of Vi with solutions decaying along V2. 

2. For < /c < 1, F = is negative, whereas A = 1 — is positive. When there is no 
damping, the point (0, 0) is a center. As fi is increased, the center transforms into a 
stable spiral for jj < 2y/l — k as shown in Figure [2|a). The frequency of spiralling is 
u ^ a/1 - A; - /i2/4. As yu 2y/l - k~ , v 0+. For /i > 2 a/1 - k, the fixed point 



(0, 0) transforms to a stable node (Figure ^b)). For /i = 2Vl - k, it is a degenerate 
node. However, degenerate nodes are borderline cases and are sensitive to nonlinear 
terms. 
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3. For /c = 1, (|5]) and ^ simplify to, 



9=6^ 

Oi = — sin 6{1 — cos 6) — fi 6i 



(9) 
(10) 



In the linearized dynamics, 6i decays exponentially as e^^*. In the phase plane, all 
trajectories move along a straight line with slope —fi and stop on reaching the 6 axis. 
However, the inclusion of nonlinear terms changes this situation. 
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(a) stable spiral, /i < 2%/! — k (b) stable node, /i > 2y/l — k 

FIG. 2. Phase trajectories around (0, 0) for k = 0.91 showing (a) stable spiral for fi = 0.3 and (b) 
stable node for = 1. 



A. Nature of (0, 0) with nonlinearities 



1. < k < 1 and /i < 2y/l - k : 

To include the effect of nonlinear terms, let us define two new variables, 

9 = rcos(f), 9i = rsm(f) (11) 
Equations ^ and ^ then may be written as, 

r = r [cos sin (p — 2y/l — k sin^ 0] — sin sin(r cos — k cos(r cos 0)] (12) 
= —Vl ~ k sin(20) — sin^ ^ sin(r cos 0)[1 — /c cos(r cos0)] (13) 



We wish to examine the fixed point (s) in the r — plane corresponding to (0, 0) in the 6 — 6 



plane, to determine their true nature. Strictly speaking, (12) and (13) are meaningful only 
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when r > 0. Neither nor (j) have any meaning when r = 0. Hence, we may assign any 



arbitrary function /(0) to 0' at r = without altering physical predictions. However, (13) 
describes accurately the approach to r = (if any) in the r — plane at arbitrarily small 
scales. We therefore set /(0) equal to the limiting value of 0' as r — )■ 0. 



f{<P) 



lim I 

r-s>0 



k k 
Vl-A;sin(20) + - cos(20) - (l - 2) 



(14) 



Equations (12) and (13) are periodic in with period vr. Hence, the phase portrait in the 
r — plane is periodic along the 0-axis with period vr. This means that in the 6 — 61 plane, 
the phase space is symmetric about (0,0). 

Using the identity (v^r^)^ + {k/2)^ = (1 - k/2Y, one can write /(0) in the form. 
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[cos(20 + a) - 1] 



(15) 



These 



where a = 2tan^^ — k, n = 0,1,2,.... Therefore, fixed points (0,0*) in the r 
plane, where /(0*) = 0, are given by 0* = nvr — tan^-*^ \/l — k with n = 0, 1, 2, . . 
correspond to the point (0, 0) in the 6 — 61 plane. 

The fixed points (0, 0*) in the r — plane are separated by nn (where n is any integer). 
Hence, in the 6 — 61 plane, there are no trajectories that can approach (0, 0) along two 
independent directions. So we can say that (0, 0) in the 9 — 9i plane cannot be a stable 
node. In the r — plane, close to some fixed point (0,0*), if there exist trajectories that 
approach this point and stop there, the corresponding fixed point (0, 0) in the 6 — 61 plane 



cannot be a stable spiral. For a spiral, — ?■ ±00 as r — )■ 0. As (12) and (13) are periodic 
in with period vr, the nature of all fixed points on the 0-axis separated by vr is identical. 
So we may choose to investigate (0,0* = — tan^-*^ \/l — k). Linearization about this point 
incorrectly predicts the whole axis to be a line of fixed points. So, we must include the 
effects of the nonlinear terms. Let e = — 0*. For small r, we may write (12) and (13) as, 

(16) 

cos(2e)] + 0(r2) (17) 



v^r^- (1 - ^)sin(2e) + 0(r 



Consider an initial condition, r(r 



0) = Tq and e(r = 0) = cq, where < eo < 



tan ^(a/I — k/2). For any finite positive value of 1 — /c, eo is finite and < 1/2. There- 
fore, ro may be chosen sufficiently small, < tq ^ eo, so as to make all terms of O(r^) and 
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higher, neghgible compared to the leading terms in (16) and (17). When these terms are 



neglected, (16) and (17) may be solved to yield, 

1 



e(r) 
r(r) 



tan 



-1 



.(2 — k)T + cot eo- 
Tq sin eo^/[{2 — k)T + cot cq]^ + 1 exp(— -\/l — k r) 



(18) 
(19) 



According to this solution, the trajectory monotonically approaches the point (0, 0) in the 
r — e plane as r — )■ oo. This behaviour will hold even with inclusion of nonlinear terms, 
provided, the terms independent of r and of 0{r) remain dominant over the entire trajectory. 
The following arguments establish that it is indeed so. 

First, the trajectory cannot reach the (p axis at a positive value of e. This is because on 
the axis, r' = and (f) = /(0) < in between two fixed points. Thus, there is already a 
trajectory running along the (p axis directed towards e = 0. 

For the initial point (ro, cq), with the choice < Tq ^ cq < 1/2, both r and e will start to 



decrease as per (16) and (17). Hence will become more negligible compared to r. Also, 



for all < e < cq, the e independent term, namely, —a/1 — kr, in (16), will be dominant. 
However, if e decays more rapidly, such that at some stage r ~ e, then the O(r^) term will 
contribute on the same scale as the first term in (17), which is O(e^). Similarly, the 0{r^) 



term will contribute on the same scale as the 2nd term in (16) if in the course of decay, at 
some point e ~ r^. However, such situations will never arise as is shown below. 

Let us assume that tq = e^^ 1 and that e decreases very rapidly, such that, at some 
instance, r = e^. As r has decreased monotonically from its initial value, we must have 



r^^"^ < Tq''^. However, r. 



1/2 




^ 1. Hence, along the entire trajectory, up to this 



instance, terms of O(r^) and higher are negligible compared to the leading terms in (16) and 



(17). Thus the solutions (18) and (19) are valid and give the correct orders of magnitudes 
of the dynamical quantities. 

As < e < ej!|° ^ 1 and < eg < 1/2, we may write, tane ~ e and tan eg ~ cq. Hence, 



< tane < tan^'^eo, which implies tane < (Vl — k/2Y^ <^ 1. Combining this with (18) 
and using the fact that (2 — /c) ~ 1, we get r > (coteo)"'^^, a very large quantity. Meanwhile 



(18) and (19) together imply, 



ro sm eo 



(20) 



-exp{\/l — k t) 

Both ro sineo and the quantity within brackets are <^ 1, which implies that r ^ e^. This is 
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in contradiction to the initial assumption that r = e^. Therefore, we conclude that starting 
with the prescribed initial condition, r will never become equal to e^. This ensures that 
along the entire trajectory, the terms independent of r and of 0(r) remain the dominant 



terms in (16) and (17). Both r and e will decrease monotonously toward their respective 



zero values. Neither can the trajectory cross the curve r = nor reach the axis before e 
becomes zero. Thus, the trajectories in the r — e plane must approach (0, 0) tangential to 
the (f) axis and slow to a halt there. 

In the 9 ~ 9i plane, the above arguments imply that, trajectories exist which start at a 



finite distance from (0, 0) and reach this point along a line of slope — -\/l — k. Also, no other 
such line with a different slope exists. These facts clearly establish that (0, 0) is a stable 
degenerate node (Figure [3]). 
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(a) stable degenerate node, /i — 0.6 (b) central region magnified 

FIG. 3. Phase trajectories around (0, 0) for k = 0.91 and /i = 2\/l — k. 

2. k = 1 : 

Proceeding as before, ^ and (10) may be written in terms of r and as, 
r = r [cos (p sin (f) — fi sin^ 0] — sin sin(r cos 0) [1 — cos(r cos 0)] 



— yU cos sm — sm — 



COS( 



sin(r cos 0) [1 — cos(r cos 0)] 



/(0) is defined as, 



/(0) = lim = —fx cos sin — sin^ 

r— s>0 



The phase portrait is periodic in with period it. The fixed points in the r — 
the form (0, 0*) are given by. 



(21) 
(22) 

(23) 
plane of 



')l = mi and 02 = nvr — tan ^/i ?t, = 0,1,2. 



(24) 
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For — tan~^ fi < (f) < 0, f{(f)) > whereas for < < tt — tan^^ fi, f{(f)) < 0. The positive 
and negative nature of /(0) repeats periodically along the (p axis. 
The Jacobian at the points (0, nvr — tan~^ yu) given by, 

J(0,^2)= I ° I (25) 

is traceless and has a negative determinant A = — yU^. So, this family of fixed points are 
saddles having stable manifold along r axis and unstable manifold along axis. 

Linear analysis of the family of fixed points (0, mr), incorrectly predicts (p = nn to be lines 
of fixed points. Let us examine the point (0, 0) in the r — (p plane for simplicity. Consider 
the condition, 

< 101 < r < mm{l,/i} (26) 
If (26) holds, then neglecting terms of O(r^) and smaller in (21) and (22), we may write. 



0r + rii (27) 

2 

~fi4>-j-^^ + r]2 (28) 



where rji ~ and 772 ~ 0^ or 0^r^, whichever is larger. Note that as long as (26) is 
satisfied, both rji and 772 can at most be of the order of r^. 

Let us take the initial condition 0o = and < ro ^ 1. Then, 0'(t = 0) = — rQ/2 and 
r' (t = 0) = 0. Hence, will start decreasing and become negative. As a result, r' will 



become negative and remain so until or r vanishes, provided (26) remains true. It is seen 
that as long as the trajectory is above the curve = — r^/2yU, (26) is satisfied and both r' 
and 0' are negative. Therefore, the trajectory approaches and eventually crosses this curve, 
where 0' is still negative, being of the order of 0^. 

Let us consider the 'trapping region' in Figure |4| which shows the phase flow on the 
curves = — r and = — r^/2/i. At any point on the line = — r for which r ^ mm{l, /i}, 

r r 

As r ^ /i, |0'/r'| ^ 1. Thus the phase flow is almost vertically upward, as shown in Figure 



4j Everywhere inside the region, (26) is satisfied and hence r < and finite. Consequently, 
after entering the region at point P, the trajectory must constantly move towards left. 
Again, it cannot penetrate the curves AP or AB, because other trajectories are actually 
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FIG. 4. Trapping Region. 

flowing inward across tliem. Hence, we liave trapped it. Upon arrival at any point on the 
arc AP, a trajectory must inevitably land up at (0,0). Note that at any point on the line 
= — mr, (0 < m < 1), 

(j)' — {m^ + l/2)r + 0{im?r'^) 

r' mr + 0{m'^r'^) 

For any finite value of m, this approaches — oo as r — )■ 0, meaning that the phase fiow is 

almost vertically upward on any line of non-zero slope near (0, 0). Therefore, the trajectory 

must reach (0, 0) along the r axis. 

Thus, the fixed points (0, mr) in the r — plane are stable nodes having slow eigenvector 

along r axis and fast eigenvector along (p axis. In between these, lie the saddle points 

(0, nvr — tan^^ /x) (Figure [s]). 




(a) stable node and saddle points (b) region near (0, 0) magnified 

FIG. 5. Phase portrait in the r — (j) plane. 



These results from the r — (j) plane mean that in the 6 — 6i plane, two trajectories exist 
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which reach (0, 0) along the hne of slope —fi and all other neighbouring trajectories reach it 
along the 6 axis. In other words, (0, 0) is a stable node here. 




(a) stable node (b) central region magnified 

FIG. 6. Phase trajectories about (0, 0) in the 6 — 9i plane for /c = 1 and fi = 0.5. 



Figure uk shows the nature of the fixed point at (0,0) over the entire parameter space. 




(a) (0, 0), curve is ii — 2^/1 — k (b) (±i7i, 0), curve is ji — l^k — \/k 

FIG. 7. Nature of fixed points (0, 0) and (iOi, 0) over the k — ^ plane 



IV. NATURE OF THE FIXED POINT (J7i,0) 

This fixed point exists when k >1. Linearization of ([s]) and ([6]) gives us the Jacobian at 
(fii, 0) as, 

J(f^i,0)= I / , ^ I (29) 
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For < < 2\/ k — J,, {ill, 0) is a stable spiral with eigenvalues given by, 



4 



Trajectories spiral in with an angular frequency v ~ — ^/k) — fi^/A, while their radial 
distance decreases as e~'^*/^. As /i — )■ 0, this decay rate vanishes and (fii,0) turns into a 



center (Figure 8). Also, u vanishes as — 2 a/A; — 1/k , representing a smooth transition 
to a stable node, similar to the behaviour of the fixed point (0,0). When fx > 2^/k — ^, 




Ol 0.0 




-0.6 -0,4 -0,2 0.0 0.2 0.4 0.6 



(a) center, fi = 



(b) stable spiral, ^ = 0.1 



FIG. 8. Phase trajectories around (0,0) for k = 1.1 and < /i < 2a/A; — I 



F^ — 4 A > 0, hence (f2i,0) is a stable node (Figure [o]), with eigenvalues and eigenvectors 
given by. 



M,2 — 7, , Vi^2 



1 



2 ■ ^ (-/x±6)/2 

where ^3 = \/ fi"^ — 4(A; — l/k). Both Ai and A2 approach the value —a/A; — 1/k as /i — )■ 



2a/A; — 1/k , indicating a stable degenerate node. For fi = 2^^ k — F^— 4A = 0. In the 
linear stability analysis, (fii,0) is a stable degenerate node with a single eigenvector. 



VA;- 1/A; 



(30) 



corresponding to the eigenvalue A = —^k — 1/k. However, degenerate nodes can be trans- 
formed into stable nodes or stable spirals due to perturbation introduced by nonlinear terms. 
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0.4 -0,2 0.0 0,2 0,4 




(a) stable node, = 1 



(b) stable node, ^ — 2.5 



FIG. 9. Phase trajectories around (0,0) for k = 1.1 and fi > 2\ k — p 



A. Nature of (Oi,0) with nonlinearities 



As discussed in subsection III A, ([S]) and (|6]) may be transformed to equations in r and 
and with the substitutions, 9 — Vti = r cos and 6i = r sin 0. Define, 



/(0) = lim 



-xlk-\ sin(20) - 1 (A; - 1 - 1) cos(20) - i(A; - 1 + 1) (31) 



/(0) is negative at all points on the axis except at the fixed points given by (0, 0*) with 



nix — a, where a = tan ^ ^Jk — 1/k, {n = 0, 1, 2, ... ), where it is zero. These fixed 

we can argue 



points are separated by nn. Then, by the same reasoning as used in III A 



that (0, 0) cannot be a stable node. The remaining possibilities are a spiral or a degenerate 
node. 

Let us consider the fixed point (0, —a), and let e = + a. Then we may expand r' and 
e' upto 0(r), 



r-' = -^(^k - ^ + 1^ { sin(2a) - sin(2e))r + 0{r'^) 

1 \ 2 3 I r ^ 

K— - + lsin e — -r\ 1 cos a 



(32) 



9 / 1 

--Wl — — cosa sin(2a) r e + O(r^) 
4 V K 



(33) 



We choose an initial point, (ro,eo), and < tq ^ cq < min{l,a}. This choice ensures 
that both r and e will start decreasing. In the course of this monotonic decay, r cannot 
reach zero before e becomes zero, as there is a straight line trajectory moving downward 
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along the axis. From (32) and (33) we note that r' and e' each contains an e independent 



term of 0{r). Therefore, as r and e decrease toward their respective zero values, the 1st 
order approximation gets even better. However, if at some stage, e ~ r, then the O(r^) 



terms would contribute on the same scale as some of the terms of 0(re) in (32) and (33). 
But the following argument rules out such a possibility. 

Let us consider a specific case and choose tq = eg". If e is to become 0(r), at some stage, 
we must have r = e^. However, for r = e^, we have 



k 



0(e' 



e = — 



1 



_\_ 

fc2 



+ 0(6^) 



whereas along the curve r = e^, dr/de = 2e. Therefore, for a given value of A; > 1, we can 
always select an eo, sufficiently small, for which, at all points on the curve r = contained 
between e = and e = eo, the ratio r'/e' > dr/de. This would guarantee that the trajectory 
cannot penetrate down this curve, which means that e cannot reach zero before r does. Thus, 
for a suitable choice of initial conditions, the trajectory must slow to a halt at (r = 0, e = 0). 
In the 6 — 01 plane, this means that there exist trajectories which start at a finite distance 
from {Qi, 0) and reach it along the line of slope —a. Also, there is no other such line with a 



different slope. Hence, (fii,0) is a stable degenerate node (Figure 10). Figure uId gives the 
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(a) stable degenerate node 



(b) region around (f2i,0) magnified 



FIG. 10. Phase trajectories for k = 1.1 and = 2\ k — l 



nature of the fixed point at (±r2i, 0) in different parts of the parameter space. 
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V. NATURE OF THE FIXED POINT (vr, 0) 

The Jacobian matrix at (vr,0) is given as, 

J(vr,0)= f ° M . (34) 

The fixed point is a saddle for all values of k and remains so even with the inclusion of 
nonlinear terms. The eigenvalues and corresponding eigenvectors are given by, 

_ -/i±6 _ ( 1 

^^1,2 — ; Vi 2 — 



-/i±6)/2 



where ^2 = V 1^^ + ^(^ + 



A. Trajectories 

Damping of the bead leads to some qualitatively different trajectories in addition to 
those observed for the frictionless cas^^. These are mainly the different kinds of damped 
oscillations (underdamped, critically damped, overdamped) about the stable equilibrium 
points. Some of these are illustrated with the following numerical plots. 




(a) k = 0.75, = 0.05, 0' (0) = (b) /c = 4, ^ = 0.5, 9' (0) = 

FIG. 11. Underdamped oscillation about (a) ^ = and (b) 9 = Qi. 



VI. PHASE PORTRAITS AND BIFURCATION 



For < A; < 1, the fixed point at (0,0) transforms its nature as the damping coefficient 
is varied. It is a center at /i = 0, as /i increases, it becomes a stable spiral. At = 2^1 — k, 
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(a) k = 0.75, n = 1.5, 6' (0) ^0 (b) A: = 4, = 4, e''(0) ^ 

FIG. 12. Overdamped oscillation about (a) 9 = (b) 9 = 



it turns into a stable degenerate node. It makes a smooth transition to a stable node as 
damping is increased further. Thus, a spiral-node bifurcation takes place at this critical 
condition (Figures |2] and [3]). Physically, as damping is gradually increased from 0, the 
system undergoes a continuous transition from undamped oscillations of the bead about 
6 = (center), to underdamped oscillations (stable spiral). At fi = 2\/l — k, the system is 
critically damped (degenerate node) and becomes overdamped (stable node) as /i is increased 
further. 




e g e 

(a) unstable spiral, = — 1 (b) center, /i = (c) stable spiral, /i = 1 

FIG. 13. Phase portraits for k = showing degenerate Hopf bifurcation. 



For negative damping, (0, 0) becomes an unstable spiral and changes to an unstable node 
as /i is made more negative. Consequently, as one crosses fi = 0, the fixed point (0,0), 
undergoes a degenerate Hopf bifurcation (Figure [l3|). 

With increase in the angular speed of the hoop (i.e., k), the stability of the origin degrades 
continuously. When = 1, (0, 0) is a weak center. A special case of Hopf bifurcation occurs. 



when /i is swept from negative to positive values acroos 0, keeping k fixed at 1 (Figure 14). 
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As k is increased beyond 1, (0,0) transforms from a stable (/x > 0) or unstable {fi < 0) 
node to a saddle. Two new stable nodes appear at ±f2i = ±cos^^(l//c) and branch out in 
opposite directions. Thus, a supercritical pitchfork bifurcation occurs at {A; = 1} (Figure 

Tsk). 
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(a) unstable node, fi = —0.2 (b) stable node, ^ = 0.2 

FIG. 14. Phase portraits for k = 1. 








=0.5 




,1 k=1.2S 



(a) Supercritical bifurcation at fc = 1, /i = (b) Supercritical bifurcation at fc = 1, /i = \/2 

FIG. 15. Section of the bifurcation diagram for (a) = and (b) /x > 0. In (a) solid curve represents 
center, dashed curve represents saddle. In (b) solid curve represents stable node, densely dashed 
curve represents saddle and sparsely dashed curve represents stable spiral. Supercritical pitchfork 
bifurcation occurs at k = 1 and spiral-node bifurcation occurs at A: = 0.5 and k = 1.28. 



As k increases from 1 to oo, Qi varies from to 7r/2. The fixed point (f2i,0), is a 
center for zero damping, a stable spiral in the region < /i < 2a/A; — 1/k (underdamped 



oscillation), and becomes a stable degenerate node at critical damping fi = l^k — Xjk. For 
the overdamped condition /i > 2^k — \/k, it is a stable node. 
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For negative damping, we get just the unstable counterparts. Accordingly, a spiral-node 
bifurcation is observed at /i = ±2 a/A; — 1/k (Figures 8, 9, 10 and 15 )b. A degenerate Hopf 
bifurcation is observed for > 1 and /i = (Figure 16). The fixed points (ivr, 0) are saddles 






(a) unstable spiral, /i = — 1 (b) center, /i = (c) stable spiral, fJ. — 1 

FIG. 16. Phase portraits for k = 4 showing degenerate Hopf bifurcation. 

for all values of k. They have saddle connections between them at /i = 0, which break in 
opposite directions for positive and negative damping. 

The above observations are summarized in Figure 17 and Table [T] below. 




FIG. 17. The bifurcation diagram 



In Fig [17} red denotes stable node, green denotes stable spiral, blue denotes unstable 
spiral, yellow denotes unstable node, brown denotes saddle, pink denotes center, gray denotes 
unstable degenerate node, and black denotes stable degenerate node. 

Up to now, we have limited attention to those bifurcations resulting from a variation of k 



or a variation of fi. From Figure 17, we see that the curves fi = 2(1 — k), /i = 2(k — l//c). 



and k = 1 divide the parameter space into 8 distinct regions of different dynamics. All these 
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TABLE I. Bifurcation Table. 



Points in parameter space 


Bifurcation along k 


Bifurcation along fi 


Figure references 


1) (fe,o) ■,k^l 




degenerate Hopf 


Figures 


13 


16 




2) {k,±2Vl -k) ■,0<k<l 


spiral-node 


spiral-node 


Figures ^ 


I 3 


15 




3) (1,0) 


supercritical pitchfork 


Hopf 


Figures 


13 






14 


, 15 


I 


4) (1,^) ; ^/O 


supercritical pitchfork 




Figures 


2 




I, 


14 3 


15 


3 


5) ik,±2^k-l) ; A;>1 


spiral-node 


spiral-node 


Figures 8 


10 


15 


3 



regions meet at the point {k = 1, fi = 0}. Traversing suitable curves in k — fi space, one 
can move from any one region to another, yielding new kinds of bifurcation. Following such 
a curve amounts to keeping a certain function a{k,fi) constant, while varying some other 
function P{k,^). Mathematically, the possibilities are rich. But whether it is possible to 
actually implement this in the bead-hoop system is subject to further inquiry. However, this 
would attain physical significance if there exists another system where a and (3 themselves 
are the control parameters. 



CONCLUDING REMARKS 



The simple introduction of damping to the bead-hoop system enriches its dynamics and 
leads to various new modes of motion and different classes of bifurcations. We have studied 
this system over the entire parameter space and presented phase portraits and trajectories. 
This serves to illustrate the qualitative changes in the system's dynamics across different bi- 
furcation curves. We have presented exact analytical treatment of the borderline cases where 
linearization fails, for which no general methods are available in the literature. The method 
of transforming to polar coordinates and using order of magnitude arguments, employed in 
this article, can serve as a useful technique for other dynamical systems as well. 
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